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FOREWORD 


This report describes a portion of the results obtained on NASA Grant 
NSG 3044, This work was done under subcontract to the University of Illinois, 
Urbana, with Prof. S.S. Wang as the Principal Investigator. The prime 
grantee was the Massachusetts Institute of Technology, with Prof. F.J. McGarry 
as the Principal Investigator and Dr. J.F. Mandell as a major participant. 

The NASA - LeRC Project Manager was Dr. C.C, Chamis. 

Efforts in this project are primarily directed towards the development 
for finite element analyses for the study of flaw growth and fracture of 
fiber composites. This report presents study of boundary layer problems in 
composites based on an exact solution. 
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ABSTRACT 


A study of boundary-layer stress singularities in multilayered fiber- 
reinforced composite laminates is presented. Based on Lekhnitskii' s stress 
potentials and the theory of anisotropic elasticity, formulation of the problem 
leads to a pair of coupled governing partial differential equations. An eigen- 
function expansion method is developed to obtain the homogeneous solution for 
the governing P.D.E's. The order or strength of boundary-layer stress singu- 
larities is determined by solving the transcendental characteristic equation 
obtained from the homogeneous solution for the problem. Numerical examples of 
the singular strength (or singular eigenvalues) of boundary-layeiT stresses are 
given for angle-ply and cross-ply composites as well as the cases of more general 
composite lamination. 
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1 . INTRODUCTION 


The response of a multilayered fiber-reinforced composite laminate near 
its geometric boundaries has been a subject of intensive investigation during 
the last decade. Both experimental studies and approximate analytical solu- 
tions have indicated that complex stress states with rapid change of 
gradients occur along the edges of composite laminates, for example, Refs [1- 
15]. This phenomenon is considered to result from the presence and inter- 
actions of geometric discontinuities of the composite and materials discon- 
tinuities through the laminate thickness. The anomaly has been found to occur 
only within very loca]. region near the geometric boundaries of a composite 
laminate, and is, therefore, frequently referred to as "boundary-layer effect" 
or "free-edge effect" — a problem unique to composite laminates and not 
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It has been shown further that 


the boundary-layer effect is three-dimensional in nature and not predictable 
by the classical laminate theory (CLT) [16,17]. The boundary-layer effect is 
apparently one of the most fundamental and Important problems in the mechanics 
and mechanical behavior of composite laminates. The high stresses developed 
in the boundary-layer region coupled with the low interlaminar strength are 
certainly of critical significance in aggravating the failure of composite 
materials and structures, For example, boundary-layer stresses have been 
observed to be responsible for the initiation and growth of local hetero- 
geneous damage in the forms of interlaminar (delamination) and intralaminar 
(transverse cracking] fracture in composite laminates under static loading 
[3,18], They are considered to have even greater effects on the long term 
strength of composite laminates under cyclic fatigue loading [19,20], 
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While the significance of boundary-layer effects has long been recognized, 
research progress on this subject has been relatively slow. The situation is 
apparently caused by the inherent complexities involved in the problem: the 

strong anisotropy of mechanical properties of each individual ply, the abrupt 
change of materials properties through the laminate thickness, the geometric 
discontinuity along laminate boundaries, and the coupling between in-plane 
and transverse deformations and stresses near the edges of the composite lami- 
nate. According to Pagano [14], analytical studies to date may be roughly 
classified into two general categories: approximate theories and numerical 

solutions. The first approximate solution for finite-width composite laminates 
was proposed by Puppo, et al. [4] based on a laminate model containing aniso- 
tropic laminae and isotropic shear layers with the interlaminar noirmal stress 
being neglected throughout the laminate. Other approximate theories were also 
attempted to examine the problem such as the extension of the higher-order 
plate theory [21] by Pagano [10], the perturbation method by Hsu, et al. [12], 
and a boundary-layer theory by Tang, et al. [11]. Recently, Pagano [14,15] 
has developed an approximate theory based on assumed in-plane stresses and 
the use of Reissner's variational principle. Even though there is no singu- 
larity involved in the formulation, the approach has certain features signifi- 
cantly Important in objectively determining detailed laminate stress fields. 

The study of edge stresses in composites by using a numerical (finite differ- 
ence) method was apparently first made by Pipes, et al.[5]. Isakson and 
Levy [6] developed a finite-element scheme containing membrane elements, 
which closely resemble the laminate model of Puppo et al. [4]. Later finite 
element studies on this subject by Wang, et al. [13] and Herakovich, et al. 

[18] led to numerical solutions similar to that given by Pipes [5]. Due to 
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the singular nature of the problem, a large number of elements, especially 
through the thickness direction, are required in conjunction with a lengthy 
extrapolation procedure in order to achieve satisfactory solutions even for 
a simple two or three layer laminate. Improved finite-element methods by 
using a more complex element Stiffness formulation baaed upon Maxwell stress 
functions [7] and by hybrid-stress elements [22] have been able to achieve an 
expedient computation with significantly less elements, Unfortunately, the 
refinements do not guarantee [23] the convergence and accuracy of the numer- 
ical solutions because of the singular nature of the boundary -layer stress 
field. That is, with each more refined analysis, numerical values of the 
maximum interlaminar stresses are shown to rise with continuously decreasing 
element size. The quest, apparently, is to show that a stress singularity 
exists at the edge of a composite laminate. 

From a linear elasticity point of view, it Is well known, that stress 
singularities are prevalent at the comers of geometric boundaries joining 
dissimilar materials (see, for example, [24-26]). Unfortunately, the search 
for the order of stress singularities for the boundary-layer region in a com- 
posite laminate containing anisotropic plies has not been successful to date, 
to the authors' knowledge. Since the singular boundary-layer stresses are 
extremely localized in nature, the precise nature of the boundary-layer effect 
will not be fully understood until the exact order of the stress singularities 
is defined. In this paper, the first in succession, a rigorous theoretical 
investigation of the free-edge stress singularity in composite laminates is 
presented . 

In the next section, a mathematical model and basic equations for each 
lamina of the composite are presented. Based on the theory of anisotropic 
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elasticity and Lekhnitskil's stress potentials [27], a pair of linear 
governing partial differential equations is derived. Appropriate 
near field boundary conditions, end loading conditions and interface contin- 
uity conditions are given also. The homogeneous solution for the problem is 
obtained in Section 3 by an eigenfunction expansion method. The solution proce- 
dure used to evaluate the exact order of the boundary-layer stress singularity 
is presented. Degenerated cases of commonly used cross-ply composite laminates 
are examined also. Numerical examples of determining the edge stress singu- 
larities for graphite-epoxy composite laminates with various fiber orientations 
are given in Section 4. As will be shown later, the existence of the free- 
edge stress singularity in composite laminates is proven mathematically in 
this paper. It settles, once and for all, the previous conjecture of boundary- 
layer stress singularities in composite materials, and provides a rigorous 
mathematical method for determining the exact value of the edge stress singu- 
larity, which is the fundamental basis for the boundary layer theory in com- 
posite materials and structures. 


2. FORMULATION 


2.1 Basic Equations 

Consider a composite laminate composed of flber-relnforced piles with 
constitutive equations described by i,i;enorallzed Hooke's law In the con- 
tracted notation as 
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(l.j - 1,2, 3, 4, 5, 6), (1) 


where the repeated subscript indicates summation and is the compliance 
tensor. The engineering strains, e^, in Eq 1 are defined In a Cartesian co- 
ordinate system by 
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where u, v and w are the components of displacements. The stresses, are 
defined in an analogous manner in the Cartesian coordinate system. 

The composite laminate considered here has a finite width and is sub- 
jected to surface tractions acting in planes normal to the generator of 
the lateral surface and not varying along the generator, i.e., along the z- 
axis (Fig. 1). The composite is assumed to be sufficiently long that, in 
the region far from the ends, the end effect is neglected by virtue of the 
Saint Venant principle. Consequently, the stresses in the laminate are in- 
dependent of the z-coordlnate. The case of a finite-width composite laminate 

subjected to a uniform axial strain, e =a, along the z-axis has been inten- 

z 

sively studied by many researchers [5-13 ] . The special case in which 
stresses and displacements are independent of z corresponds to the well known 


5 



generalized plane deformation [27] . Under these assumptions, the equations of 
equilibrium without body ;. ;fce read 
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Following the procedure in [27], it can be shown after some mathematical 

manipulation that the general expressions for displacements and the stress 

component a have the following form: 
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1 2 
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(4a) 


(4b) 


(4c) 

,4,5,6). 

(4d) 


The unknown functions, U, V and W, depend on x and y only, and can be shown 
easily to obey the following relations : 




ft V - 

W Vj +S23(A3X + A3y + A3), 

If “ ®53'V V + V- 

3y ■ ®6j“l ^ 2 '' ^ ^3^’ 


where 


^ij “ ®ij " ^13®j3'^^33’ 


(i, j = 1 , 2 , A, 5, 6 ) . 


(5a) 

(5b) 

(5c) 

(5d) 

(5e) 

(5f) 


It is obvious that the constants, u , v , w and w (i = 1,2,3) in Eqs 4a-4d 

000 1 

characterize the rigid-body translations and rotations of the solid. A^^ and 
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represent the bending of the laminate in the x ~2 and y-z planes. 
characterizes the uniform axial extension of the composite laminate, and 
A^, the relative angle of rotation about the z-axia. 

2.2 Governing Partial Differential Equations 

Introducing Lekhnitskii’s stress potentials, F(x,y) and fCxjy) [271, 
such that 
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one can show that the e luations of equilibrium are satisfied identically. 
Eliminating U and V from Eqs 5a, 5b and 5e and W from Eqs 5c and 5d by dif- 
ferentiation, we obtain the following system of governing partial differential 
.aations for the problem; 
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I L^F + 1. f » 0, (7b) 


where L 2 > and are line'*'., differential operators defined as 
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2.3 Boundary and End Conditions 

Assuming that the edges of a composite laminate, 3 Bp, are traction free 
and the interface of the mth and Cm+l)th plies is a straight line meeting 
the traction-free edge at a right angle (Fig. 1), one can obtain the following 
boundary conditions along 3 Bp: 


0 ■ x a •X 

X xy xz 
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( 8 ) 


The conditions at the ends of the composite laminate may have the foinn from 
the statically equivalent loads as 


• 0 0 
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where the integrals are taken over the entire area B of the cross section, 

and P , M , M and M are the applied force, bending momenta and twisting 
z X y t 

moment acting on the ends, respectively. 


2. A Interface Continuity Conditions 

Consider a portion of the laminate cross section composed of the wth and 
(m+l)th fiber-reinforced laminae, as shown in Fig. 1. Assuming that the 
plies are perfectly bonded along the interface 3Bp, one can immediately es- 
tablish the continuity conditions of the stresses and displacements along 
the interface as the following: 
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and 
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where the superaerlpta denote the mth and (w-»*l)th plies in a composite 
laminate, and i\^ and n^ are components of unit outward normal to the inter- 
face. 
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3 . HOMOGENEOUS SOLUTION AND FREE-EDGE STRESS SINGULARITY 


The governing equations, 7 a and 7 b, are coupled, linear partial differen- 
tial equations with constant coefficients related to the anisotropic elastic 
constants of each individual lamina. With the aid of aforementioned noar-fiold 
boundary conditions and interface continuity conditions, the homogeneous solu- 
tion for the governing P.D.E.'s can be determined easily. The homogeneous 
boundary conditions and interface continuity conditions also provide the in- 
formation for determining the important strength or order of the free-edge 
stress singularity in a composite laminate, which is the major concern in 
this paper. 

According to Lekhnitskii [ 27 ], the homogeneous solution for the governing 
partial differential equations has the general form as 

F(x,y) » I F (x + u.y), H'(x,y) " I Hj. F^^(x + u^^y), (lla-b) 

k-1 k“l 


where the prime (') in Eq 11 b denotes differentiation of the function 

F^^(x ^ with respect to its argument, and the coefficients are the roots 

of the following algebraic characteristic equation 
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It can bo shown that Eq 12a cannot have a real root (thus, have to appear as 
complex conjugates), and are analytic functions of the complex variables 

•• X + Mj^y « r(e^® + Xj^e"^®)/(1 + Xj^) with X^ - (1 + iu^^)/(l - iUj^) and 
r and 0 being components of polar coordinates. Substituting the expressions 
of F(x,y) and 'l'(x,y), Eqs 11a and 11b, into Eqs 6a-6e, the homogeneous solu- 
tions for stresses may be expressed in terms of 
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where 
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We now choose Che form of as 
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Eq 15 into Eqs 13 and 14 gives 
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where the overbar denotes the complex conjugate of the associate variable. 

For convenience, we drop the superscript h associated with the above homogeneous 
solutions for stresses and displacements in this paper. 

The homogeneous solutions are required to satisfy the homogeneous boundary 
conditions and interface continuity conditions. This leads to a standard 
eigenvalue problem for determining the values of 6 . It is noted that 5 general- 
ly appears as a set of complex conjugates, which enable to make Eqs 16 and 17 
real functions by superposition. Furthermore, the value of 5 is required to 
satisfy the condition 

Re[6] > -1 (18) 


to ensure the finiteness of displacement components at the origin, where Re 
represents the real part of 6 . 

To expedite further developments, we transform the stress and displace- 
ment components from Cartesian coordinates to polar coordinates. Thus, we 
have 
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where Zj^ are defined in the polar coordinates and (j » 1,2,... 8) are 
functions of ri|^, Pj^, Pj^, t^ and 0 given in Appendix 1. 

The traction-free boundary conditions, Eqs 8a-8c, along the free- Iges 
of thewth and (m+l)th plies in polar coordinates read 
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More explicitly, the homogeneous boundary conditions, Eqs 21a and 21b and the 
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Substituting Eq 25 into Eq 22a gives 
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Equations 22b and 26 constitute a system of homogeneous linear algebraic 
equations in . The existence of a nontrival solution for 

quires the vanishing of the coefficient determinant 


re- 


|4(6)| = 0, 


(27) 


where A (6) is a six by six matrix involving 6 in a transcendental form. Thus, 
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Eq 27 is a transcendental characteristic equation for the standard eigenvalue 
problem. It has a very complicated structure as can be seen from the coeffi- 
cients of ^ in Eqs 22 and 26, and the detailed expression for A (6) is not 

given here. The investigation of the characteristic equation requires the 
employment of standard numerical techniques such as the Muller’s method [28] 
with the aid of a digital computer. The eigenvalues 6^ obtained from the 
numerical solution of Eq 27 give important information concerning the behavior 
of the edge stresses and displacements. Due to the positive definiteness of 
strain energy of the elastic body and the condition in Eq. 18, the eigen- 
value of 6 bounded by 


-1 < Re[6 ] < 0 
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(28) 


characterizes the order of the inherent singularity of the boundary- layer or 
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4. DEGENERATED CASES - CROSS-PLY COMPOSITE LAMINATES 


In Che case of a crosa-ply composite laminate, i.e., laminae with 0° and 
90® fiber orientations only, the in-plane stresses and the out-of-the-plane 
shear stresses are uncoupled by virtue of the material symmetry in each lamina. 
We shall concentrate our study on the four stress components of interest here, 
Uy, o^. The general expressions for displacements and may be 

simplified as 
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Following the same procedure shown in the previous section, the govern- 
ing partial differential equations are uncoupled and may be written as 
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The homogeneous solution for Eq 30a may be obtained in a simplier form, 
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F(x,y) » J F^(x + 
k^l 


where u, are the roots of the following algebraic equation: 

cC 


£,,(y) = + C2S^2 ^22 “ 


'12 ■ “66' 


(32) 
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Tile homogeneous stress and displacement solutions are then given as 


*4 H 


k"l 


k»l 


xy 

U(x,y) 

where 


-I F|;(x + u y) 
k«l 


"T 

I Pk ^k^^* 

k-1 ^ 


r ' 

v(x,y) " I qj^ + Mj^y) , 


k-1 


^k " ^ll^k ^12’ 


\ *• ^2l\ 222/^k- 


(33) 

(34) 


(35) 


We shall choose the form of F^(x,y) as 

F^(Z^) - C^Z^V[(6+l)(S+2)). (36) 

where C, and 3 are, as before, arbicrary complex constants to be determined 
ic 

later. Imposing the homogeneous traction-free boundary conditions along the 
free edges and the continuity condition along the ply interface, one can 
precede with the same procedure outlined in the previous section. Then, the 
eigenvalues and eigenfunctions can be determined in a manner similar to those 
in the previous cases. 
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5. NUMERICAL EXAMPLES 


From the structure of the governing p.. cial differential equations and 
the homogeneous solution for the problem, it is clearly seen that the asymp- 
totic stress and strain fields in the vicinity of the edge are governed by 
the singular terms with the strength of stress singularity 5^ determined 
from the eigenvalue analysis. Exan. ning the structure of Eqs 22b and 26, it 
is obvious that the eigenvalue solutions, and therefore, the edge stress 
singularities, are related to laminar constitutive properties and fiber ori- 
entations of adjacent plies. 

Consider a composite laminate with ply properties typical of those used in 

earlier studies [5] (values of ^Lz“ .85 Vlt found to 

make only a few % difference in 6n) : 


- 20 X 10^ psi, 


- E “ 2.1 X 10^ psi, 
T z 


G =G =G =0.85x inf' . 
LT Lz Tz 


\) MM a\) »0. 21 • 

LT Tz Lz 


(37) 


where the subscripts, L, T, and z refer to the fiber, transverse, and thick- 
ness directions of an individual ply, respectively. The influence of material 
properties of composite plies on the boundary layer stresses may be related to 
the roots of the characteristic equation, Eq 12a. With the lamina proper- 
ties given above, the roots of the characteristic equation for the graphite- 
epoxy lamina of different fiber orientations 0 are shown in Table 1. It 
appears that all the six roots y^ are purely imaginary by virtue of the 
material properties of Eqs 37. Furthermore, the y^ for +Q ply is the 


same as those for due to the in-plane rotation of fiber directions. 


Based on the material constants, y^, p^, q^^ and t^ obtained for the 
graphite-epoxy, the transcendental characteristic equation, Eq 27, can 
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Table 1 


Roots Uj^ of Characteristic Equations for Graphite-Epoxy 
Cowposite System with Fiber Orientation 0 


±0 

*^1,2 


JHa 


*"5^6 

15“ 

±0.88782 

i 

±1.10301 

i 

±1.56776 

30“ 

±0.86222 

i 

±1.06956 

i 

±2.53630 

45“ 

±0.80902 

i 

±1.03503 

i 

±3.44438 

60“ 

±0.73316 

1 

±1.01323 

i 

±4.15870 

75“ 

±0.66324 

i 

±1.00294 

i 

±4.61201 


*0is the angle measured counterclockwise from the positive 
z-axis to the fiber direction 
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be oolved numerically to provide the eigenvalueo for the horaogeneouG 
solution. For illustrative purposes » the first twelve non-integcr eigen- 
values associated with the stress solution for the free edge of a 45®/-45° 
graphite-epoxy are shown in Table 2. Eigenvalues 6^ smaller than -1 are ex- 
cluded for the reasons given in the previous section. It is seen that there 
exists one and only one eigenvalue which satisfies the required constraint 
condition, Eq 28, for this case, i.e., 6 ^ ■ -Q. 02557. In fact, the eigenvalue 
is the strength or order of the free-edge or boundary-layer sti*ess singu- 
larity, which is of major concern in this study. The fact that theie is 
only one which meets Eq 28 is observed in all cases with various other fiber 
orientations studied; the only difference is that, for each case, 6^^ possesses 
a different value. Higher-order eigenvalues occurring as integers (including 
zero) and as complex conjugates always exist and should be included for de- 
termining the complete solution when remote boundary conditions are matched 
by a numerical method to be discussed in an associated report [29], 

For the commonly used (±0) angle-ply graphite-epoxy composite, as anti- 
cipated, the order of the bou 'dary-layer stress singularity is a function of 
the fiber orientation 0. Numerical tesults of 6^ for each of the 1*5 fiber 
orientations are calculat d and shown in a graphic form in Fig. 2. It is 
clearly seen from the Figure that the composite free edge associated with the 
laminate of an approximately (±51*) ply orientation possesses the strongest 
boundary-layer stress singularity. As the 0 changes to either directions, 
the order of the stress singularity 6^ decreases rapidly. Its value con- 
verges to zero for the cases of 0 ■ 0* or 90*, since the two plies become 
identical with orthotropic elastic properties. 
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Table 2 


First Twelve Non-Integer Eigenvalues* for Free-Edge Stress 
Solutions in (±45®) Graphite-Epoxy Composite 


-2.5575658 E-2 
8.8147184 E-1 ± 
1.5115263 E 0 ± 
2.3389433 E 0 ± 
3.0913532 E 0 ± 
3.9520023 E 0 ± 
4.7440929 E 0 ± 
5,6021457 E 0 ± 
6.3962635 E 0 ± 
7.2565174 E 0 ± 
8.0497237 E 0 ± 
8.9120567 E 0 ± 


i 2.3400497 E-1 
i 7.9281732 E-1 
i 1.1158402 E 0 
i 1.7360464 E 0 
i 2.0287146 E 0 
i 2,5683871 E 0 
i 2.85SS510 E 0 
i 3.3652707 E 0 
i 3.6575937 E 0 
i 4.1479983 E 0 
i 4.4407609 E 0 


* 


Integers* 0,1,2 


.n, are always eigenvalues obtained from Eq 27 
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Fig. 2. Strength of Boundary-Layer Stress Singularity in 
(0/-0/-Q/0) Graphite-Epoxy Composites 






In the caae of a froe edge associated with two plies of more general fiber 
orientations, instead of the symmetric -HD/-© configuration, solutions for the 
eigenvalues are obtained also. To Illustrate the nature of the eigenvalues 
for this caae, wo examine the froe edge associated with 30°/© fiber orienta- 
tions in a graphite-epoxy composite, where © varies from 7.5° to 82.5°. The 
first few non-integer eigenvalues for various 0's in these cases are given in 
Table 3. The integers (including zero) are also eigenvalues, but not included 
in the Table. The case of 30°/30° graphite-epoxy is not included either since 
the two plies are identical. Again, it is observed from the Table that there 
exists only one 6^ which meets the requirement of Eq 28 and gives the dominant 
singular stress state at the edge of each of the 30°/© composite laminate. 

The degenerated case of cross-ply composite laminates which are discussed 
in the previous section is investigated also. The eigenvalues for 
the fres-edgs stresses in a graphite-epoxy composite with 0°/90“ lamination 
are given in Table 4. The dominant stress singularity in the present case has 
an order of magnitude similar to those in (±0) angle-ply composites and in more 
general (0j^/02) laminates. It is noted that the orders of the boundary-layer 
stress singularity for both angle-ply and cross-ply composites are generally 
much weaker than those associated with other typical elastostatic singular 
stress problems such as elastic crack problems. The relatively weak singularity 
for the boundary-layer stresses introduces some unique features as well as 
difficulties for the evaluation of the boundary-layer effects in composite 
laminates, which will be discussed in an associated report [29]. 
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Table 3 


"fi -"rt 

First Five Non-Integer Eigenvalues for Free-Edge Stresses 
Associated with (30° /0) Graphite-Epoxy Composite 


0 

«3. 

62 

63 

64 

«5 

7.5° 

-2.792991 E-3 

9.887419 E-1 

1.930568 E 0 

3.095930 E 0 

3.862698 E 0 



±i 8.435637 E-2 

±i 3.319540 E-1 

±i 8.103959 E-2 

±i 7.204410 E-1 

15° 

-1.986115 E-3 

9.943591 E-1 

1.949288 E 0 

2.726071 E 0 

3.441129 E 0 



±i 6.113543 E-2 

±i 2.882098 E-1 


±i 5.183480 E-1 

22.5° 

-6.725306 E-4 

9.989196 E-1 

1.989893 E 0 

2.649382 E 0 

S. 342300 E 0 



±i 2.920556 E-2 

±i 1.661871 E-1 


±i 5.276117 E-1 

30° 



— 



37.5° 

-9.077774 E-4 

1.001306 E 0 

1.984327 E 0 

2.182540 E 0 

3.204760 



±i 5.180632 E-3 


±i 1.910473 E-1 

±i 8.009168 E-1 

45° 

-3.771122 E-3 

9.624687 E-1 

1.061916 E 0 

1.821684 E 0 

2.021985 






±i 4.926242 E-1 

52.5° 

-8.292252 E-3 

9.291682 E-1 

1.168571 E 0 

1.621302 E 0 

1.906136 E 0 






±i 6.212519 E-1 

60° 

-1.367358 E-2 

9.010723 E-1 

1.354059 E 0 

1.825695 E 0 

3.518516 E 0 




±i 1.864060 E-1 

±i 6.942274 E-1 

±i 1.413627 E 0 

67.5° 

-1.897304 E-2 

8.810614 E-1 

1.321793 E 0 

1.773577 E 0 

3.488780 E 0 




±i 3.224246 E-1 

±i 7.358765 E-1 

±i 1.428440 E 0 

75° 

-2.343868 E-2 

8.719040 E-1 

1.303673 E 0 

1.743059 E 0 

3.485411 E 0 




±i 3.820090 E-1 

±i 7.550141 E-1 

±i 1.425960 E 0 

82.5° 

-2.668523 E-2 

8.778236 E-1 

1.310388 E 0 

1.725286 E 0 

3.513604 E 0 




±i 3.866969 E-1 

±i 7.501905 E-1 

±1 1.387356 E 0 


Integers, 0,1,2 n, are also eigenvalues. 




'''•M !S 
QUALITY 


Table 4 

First Twelve Non-Integer Eigenvalues for Free-Edge Stresses 
in Cross-Ply Graphite-Epoxy Composite* 


-3.33888 

E- 

2 





8.80268 

E- 

•1 





1.41674 

E 

0 

+ 

i 

3.93303 

E-1 

1.65345 

E 

0 

+ 

i 

6.85523 

E-1 

2.83449 

E 

0 

+ 

i 

1.76219 

E 0 

3.75294 

E 

0 

+ 

i 

1.1853E 

E 0 

4.29235 

E 

0 

+ 

i 

2.66884 

E 0 

5.70726 

E 

0 


i 

3.57190 

E 0 

5.79010 

E 

0 

+ 

i 

1.52461 

E 0 

7.12293 

E 

0 

+ 

i 

4.48145 

E 0 

7.81068 

E 

0 

+ 

i 

1.76401 

E 0 


^Integers, 0,1, 2, 3..., are also eigenvalues 
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6. SUMMARY AND CONCLUSIONS 


A stvidy of boundary- layer stress singularity in both angle-ply and cross- 
ply composite laminates has been presented. Formulation of the problem is 
based on Lekhnitskii’ s comp lex- variable stress functions and basic relation- 
ships in the anisotropic elasticity theory. An eigenfunction expansion method 
has been developed to obtain the homogeneous solution for the coupled govern- 
ing partial differential equations for the problem. Angle-ply and cross- 
ply composites as well as more general laminates have been studied. Tl\e 
strength of boundary-layer stress singularity for each case has been determined 
to illustrate the fundamental nature of the edge effects in composite materials. 

Based on the information obtained, the following conclusions may be drawn: 

1. Boundary- layer or free-edge stress field in a composite laminate 
is inherently singular in nature due to the geometric and material 
discontinuities . 

2. The order of boundary-layer stress singularity can be determined 
by solving for the transcendental characteristic equation ob- 
tained from the homogeneous solution of the governing partial 
differential equations. 

3. The boundary-layer stress singularity depends only upon material's 
elastic constants and fiber orientations of adjacent plies in 
composite laminates. 

4. For angle-ply and cross-ply composites as well as more general 
laminates the order of boundary-layer stress singularity is very 
weak in general. In a graphite-epoxy system, for example, 5^^ is 
much smaller than other kind of singular stress problems in 
elastostatics such as elastic crack problems. 
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APPENDIX 1 


Expressions for Equations 19 and 20 


^Ik " (Uj^Bin9 + cos9)^ 



